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Interactions between nanoscale semiconductor structures form the basis for 
charge detectors in the solid state. Recent experimental advances have demon- 
strated the on-chip detection of single electron transport 18 through a quantum 
dot (QD). The discreteness of charge in units of e leads to intrinsic fluctua- 
tions in the electrical current, known as shot noise 9 . To measure these single- 
electron fluctuations a nearby coherent conductor, called a quantum point con- 
tact (QPC), interacts with the QD and acts as a detector 4 8 . An important 
property of the QPC charge detector is noninvasiveness: the system physically 
affects the detector, not visa-versa. Here we predict that even for ideal nonin- 
vasive detectors such as the QPC, when a particular detector result is observed, 
the system suffers an informational backaction, radically altering the statistics 
of transport through the QD as compared to the unconditional shot noise. We 
develop a theoretical model to make predictions about the joint current probabil- 
ity distributions and conditional transport statistics. The experimental findings 
reported here demonstrate the reality of informational backaction in nanoscale 
systems as well as a variety of new effects, such as conditional noise enhance- 
ment, which are in essentially perfect agreement with our model calculations. 
This type of switching telegraph process occurs abundantly in nature, indicating 
that these results are applicable to a wide variety of systems. 

Noise is generally due to randomness, which can be classical or quantum in nature. 
Telegraph noise, where there is random switching between two stable states 10 , originates 
from such diverse origins as thermal activation of an unstable impurity 11-13 , nonequilibrium 
activation of a bistable system 14-16 , switching of magnetic domain orientation 17-19 , or a 
reversible chemical reaction in a biological ion channel 20 . 

In nanoscale conductors, where charge motion is quantum-coherent over distances com- 
parable to the system size, shot noise and telegraph noise have recently been shown to be 
two sides of the same coin 6 ' 7 ' 21,22 . A QD is sufficiently small that it is effectively zero- 
dimensional, and behaves as an artificial atom, holding a small number of electrons. Figure 
1(a) shows the sample used in the experiment reported here. The QD is marked by the 
dotted circle 23 . An extra electron can tunnel into the QD from the source lead (S), stay 
in the QD for a random amount of time, and then tunnel out into the drain lead (D) if 
the applied voltage bias exceeds the temperature. This single-electron transport produces a 
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fluctuating electrical current. In order to detect the statistical properties of this current, a 
sensitive electrometer with a bandwidth much higher than the tunneling rates is required. 
The electrometer is a nearby QPC that is capacitively coupled to the QD via the Coulomb 
interaction. The voltage biased QPC detector transports many electrons through a narrow 
constriction in the surrounding two-dimensional electron gas (represented with an arrow). 
The resistance of the QPC is susceptible to changes in the surrounding electrostatic envi- 
ronment, and can therefore be used to sense the presence (or absence) of an extra electron 
on the QD 4 . When the extra electron tunnels into or out of the QD, the current I flowing 
through the QPC switches between two different values [see Fig. 1(b)]. Therefore, the shot 
noise of the QD current J (randomness in the number of switches in a given time interval) is 
intimately linked with telegraph noise in the QPC current / (randomness of duration time 
in each current value). 

We now discuss the experimental procedure. From a measurement such as shown in 
Fig. 1(b), one can directly determine the rates for electrons tunneling into and out of 
the QD 6,25 . The tunneling rates are controlled by tuning the voltages applied to gates Gl 
and G2. The data presented here was taken at two different gate voltage configurations. 
Configuration A is characterized by Ti = 160 Hz, T 2 = 586 Hz, and configuration B is 
characterized by Ti = 512 Hz, T 2 = 345 Hz. For each configuration, we collected traces 
of length T = 700 s, containing around 10 5 tunneling events. More details about the 
sample characteristics and the experimental methods can be found in the Supplementary 
Information. 

In the following, we develop a theoretical model for this experiment. Ideal switching 
of the measured detector signal between two noiseless values Ii^ are identified with the 
two current levels experienced by the QPC. When an electron enters the QD, the current 
switches from Ii to 7 2 , and when the electron leaves the QD, the reverse switch happens. 
The number of "down" or "up" switches M in a given time trace of duration £, is identified 
with the number of different transport electrons that occupy the QD in that time interval, 
naturally defining a QD current variable J = M/t (we set e = 1 to count in single electron 
charge units). The analogous number of electrons TV passed by the QPC in this same time 
interval defines the QPC current variable / = N/t. The assumption of noiseless current levels 
implies that Ii < I(t) < 7 2 , while the unidirectional nature of the QD transport implies that 
< J{t) < oc (see Supplementary Information for justification). Stochastic, statistically 
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FIG. 1: Nanodevice and current fluctuations, (a) Quantum dot/quantum point contact 
structure fabricated using scanning probe lithography 23 on a GaAs/Alo.3Gao.7As heterostructure 
containing a two-dimensional electron gas 34 nm below the surface. Ballistic electron flow through 
the QPC is indicated with an arrow, and the quantum dot is indicated with a dotted circle. 
Individual electron current is induced by applying a voltage bias between source (S) and drain (D). 
(b) The flow of single electron current is deduced by the presence of switching in the QPC current, 
which acts as an on-chip electrometer, (c) Probability distribution of charge transmitted through 
the QPC, plotted on a logarithmic scale. Data set A is represented by blue asterisks and data set 
B is represented by open red boxes. The solid lines are theoretical prediction of a universal ellipse 
for both configurations, with no fitting parameter. Slight deviations from the ellipse at the ends 
of the interval are due to small-amplitude noise fluctuations in the individual current levels 24 . 

independent quantum tunneling into and out of the QD is described with rates r lj2 21 . For 
later convenience, we define the average and difference variables I Q = (h + h)/^^o = 
(Ti + r 2 ) /2, AI = ( J 2 - 7i)/2, Ar = (r 2 - r^/2. This model is capable of describing a host 
of phenomena in many fields of science. 

Taken alone, each side of the random process may be characterized mathematically with 
the probability distributions P(I : i), P(J, t) of finding a given number of electrons transmit- 
ted in a given time, or equivalently, all current cumulants ((J 171 )) (see Supplementary 
Information for a discussion of these statistical quantities). This catalog of cumulants gives 
a unique signature of the particular electronic conductor, and is also known as full counting 
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statistics 26 . For example, the first two cumulants are the average current (J), and the shot 
noise power ((I 2 )) = J dt(5I(t)5I(0)), where 51 = I(t) - (I). For the QPC and QD, the 
average current (/), (J) and shot noise power ((/ 2 )), ((J 2 )), are given respectively by 

(/) = (hr 2 + /aro/pro), «/ 2 » = (A/^rxiyrg, 

(j) = r 1 r 2 /(2r ), ((j 2 )) = <j>(r 2 + r 2 )/(4r 2 ). (i) 

In the limit of small switching rates Ti 2 —> the current and noise of the QD vanish, while 
the noise of the QPC actually diverges. 

However, a simple specification of the counting statistics of each conductor individually 
misses the important fact that the two conductors are strongly correlated by the Coulomb 
interaction between them. The simplest measure of the correlation between the systems is 
the cross-correlation ((U)) = / dt(5I(t)5J(0)), given by 

((H)) = (A/Ar)r 1 r 2 /r 2 . (2) 

This correlator is approximately constant under scaling of the switching rates, compromising 
between the behavior of either noise individually (1). Result (2) implies that the two currents 
may be either positively or negatively correlated, depending on the sign of AVAL This effect 
has a simple physical interpretation: Taking AI > 0, if T 2 > Ti, then the system typically 
spends more time in state 1 than in state 2. The current J is increased by adding another 
"up" and "down" switch to the current trace. This new event typically divides a long time 
interval spent in Ii into two segments, subtracting a short interval from Ii and adding it to 
7 2 , thereby increasing the current /. The same argument applied with T 2 < Ti leads to a 
decrease in I given an increase in J, explaining the sign of the cross-correlation function. 

Going beyond the first two cumulants, the full QPC current distribution was predicted 
to have an elliptical shape 24 as a function of the current 

\ogP(i)/t = -(g l -g 2 ) 2 /(2Ai), (3) 

where Qi^ = y / T^|7 — 7 2 ,i|- This prediction is experimentally confirmed in Fig. 1(c) for 
data sets A (blue asterisks) and B (red boxes). We stress that the solid ellipses are taken 
directly from Eq. 3 with no fitting. The generating function for the QD statistics was also 
found 27 , and the first few QD current cumulants were experimentally verified 21 . 

In order to specify the statistical correlation between the conductors, we introduce the 
joint counting statistics of both conductors. More specifically, the correlations may be 
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FIG. 2: Joint probability distributions. The logarithm of the joint probability distribution of 
detecting QPC current (x-axis) and QD current (y-axis) is given as a color density plot, where red 
indicates high probability and blue indicates low probability (a), (c) Experimental construction 
for data set A and B respectively The experimental probability distributions were generated by 
splitting the data into a large number of subtraces, each containing on average seven tunneling 
events, (b), (d) Theoretical prediction for data set A and B respectively see Supplementary 
Information for an explicit formula. 



quantified by the joint probability distribution P(7, J, t) of finding current I and current 
J in a time t (equivalently, all cross cumulants {{I n J m ))) ) or may also be specified by the 
conditional distribution functions P{I\J) or P( J\I), the probability of observing one current, 
given an observation of the other. These distributions are all related to one another by 
P(J, J) = P(I\J)P(J) = where the last equality is an expression of Bayes' 

theorem. The conditional distributions (and their associated conditional cumulants) play a 
key role in the informational approach to detection 28 . 

From the model described above, the complete characterization of the system/detector 
fluctuation statistics may be obtained from conditional master equation formalism. The 
statistical cumulants of the current fluctuations are given with the help of a generating 
function if(A, %), such that the cross-cumulants are given simply by taking derivatives, 
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((I n J m )) = d%d™H(0,0). From the mathematical derivation given in the Supplementary 
Information, we find that the joint generating function is given by 



if (A, X ) = Afo - r + \] (AAf - Ar) 2 + I\r 2 exp x, 



(4) 



generalizing previous results 24 ' 27 . The function i7(0,x) generates the current cumulants of 
the QD 27 , while H(X, 0) generates the current cumulants of the QPC 24 . Results (1) and (2) 
follow from (4). The joint generating function (4) is directly related to the joint probability 
distribution of measuring current / and current J [see Supplementary Eq. (6)]. The loga- 
rithm of this distribution has been measured and given in Fig. 2(a), (c) for configurations 
A and B respectively. The theoretical prediction for this quantity is given in Fig. 2(b), (d) 
with striking agreement. 

Having described the joint statistical properties of both currents, we now return to the 
detection question. It is important to distinguish between physical backaction from statis- 
tical/informational backaction. The noninvasive QPC detector changes its physical current 
state depending on whether or not the QD is occupied by an extra electron, while the physi- 
cal dynamics of the QD is unaffected by the state of QPC. However, by observing a particular 
outcome of the detector variable (/), this leads to conditional (Bayesian) backaction of the 
detector on the system variable (J). This informational backaction, or constrained random- 
ness, manifests itself in a variety of novel effects. We introduce the concept of conditional 
counting statistics: The statistical current fluctuations of one system, given the observation 
of a given current in the other. These statistics may be calculated from the joint generat- 
ing function (4) (see Supplementary Information), giving a mixed generator Hi(I,x) °f the 
(normalized) conditional statistics of the QD, given the observation of a current /, 



Taking derivatives, all conditional cumulants are given by ((J n )) c = ft/2 n , yielding a set of 
universal semi-circles as a function of the current /. At the endpoints of the interval, the 
QPC current is observed to remain in Ii or I 2) and therefore there can be no QD current, 
or any associated noise. The conditional QD current cumulants all have a maximum at 7 - 
The conditional current maximum ^/T\T\/2 is always larger than the unconditional current 
(J). Fig. 3 (a), (b) compare the experimental values of the first two conditional cumu- 
lants to the theoretical semicircles, with excellent agreement. The distribution described by 
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the conditional cumulants (5) is a Poissonian distribution with a generalized rate and 
effective charge e* = e/2, showing a radical change when compared with the unconditional 
distribution 27 . 

Turning our perspective around, we can pose similar questions about the conditional de- 
tector statistics, given an observation of the system current J. While the generating function 

(A, J) for these statistics may be found analytically (see Supplementary Information), we 
focus on the conditional current (7) c , and the conditional noise ((/ 2 )) c , given by 

where S = \J J 2 + (Ar) 2 . These conditional cumulants are experimentally calculated and 
compared with Eq. (6) in Fig. 3(c), (d). As a function of the QD current J, the conditional 
current tends to either Ii or I 2 as J —> 0, depending on the sign of Ar. This corresponds 
to the most likely detector current configuration in the event of no switches observed: the 
QPC current stays in one value, also implying that the system becomes noiseless in this 
limit. This is easily seen in (6) because ((I 2 )) c is proportional to J. The exception to 
this rule is the perfectly symmetric situation I\ = T 2 , where the QPC conditional average 
current is Iq. This situation corresponds to rare symmetric switching between the states, 
whose effective rate is the conditional QD current J. The corresponding QPC conditional 
noise actually diverges in this limit, because the effective switching rate is vanishing. This 
effect, where the noise in one system (monitored by another) can be dramatically larger 
than the unmonitored noise, we refer to as conditional noise enhancement The same effect 
persists in the asymmetric situation, and the maximum of the conditional noise occurs at 
J 2 = (Ar) 2 (v / 5 — l)/2. In order that the conditional noise peak exceed the unconditional 
noise, the ratio R = |Ar|Fir 2 /rjj must be less than [(>/5 - l)/2] 5 / 2 w 0.3. For data set A 
and B, Ra ~ 0.38 and Rb ~ 0.19, so only data set B exhibits conditional noise enhancement. 
In the opposite limit, J —> oc, the conditional current tends to I , and the noise tends to 
zero. This situation corresponds to rapid symmetric switching between the current states, 
whose effective rate is again controlled by J. In both limits, the typical dynamics of the 
telegraph process gets completely taken over by the transport condition. 
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FIG. 3: Conditional current cumulants. (a) Conditional QD current ( J) c plotted as a function 
of /. (b) Conditional QD noise 

((J 2 ))c plotted as a function of /. Solid line is theory from Eq.(5) 
describing universal semi-circles. Solid blue dots and open red squares denote data set A and B 
respectively, (c) Conditional QPC current (I) c plotted as a function of J. (d) Conditional QPC 
current {{I 2 )) c plotted as a function of J. Solid line in (c) and (d) is theory, given by Eq. 6. 
In (d) the horizontal dashed lines indicate the unconditional noise level, illustrating the effect of 
"conditional noise enhancement" for data set B described in the text. The experimental statistics 
were generated by splitting the data into a large number of subtraces, each containing on average 
four tunneling events. 
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Supplementary Information 

"Conditional statistics of electron transport in interacting conductors" 

E. V. Sukhorukov, A. N. Jordan, S. Gustavsson, R. Leturcq, T. Ihn & K. Ensslin 
Sample Fabrication 

Figure 1(a) shows the sample used in the experiments reported here. The structure was 
fabricated using scanning scanning probe lithography on a GaAs/Al .3Ga .7As heterostruc- 
ture containing a two-dimensional electron gas 34 nm below the surface. 

Sample Characteristics 

The measurements were performed in a 3 He/ 4 He dilution refrigerator with an electron 
temperature of about 190 mK, as determined from the width of thermally broadened 
Coulomb blockade resonances. The charging energy of the QD is 2.1 meV and the mean 
level spacing is 200 — 300 fieV. 

Experimental Method 

The conductance of the QPC, Gqpc ) was tuned close to 0.5 x 2e 2 /h. We apply a dc bias 
voltage between the source and the drain of the QPC, Vqpc = 250 \iV and measure the 
current through the QPC, Iqpc, which depends on the number of electrons TV in the QD. 
The current signal was digitized with a sampling frequency of 100 kHz, thereafter software 
filtered at 4 kHz using a 8th order Butterworth filter. To avoid tunneling due to thermal 
fluctuations, the QD was voltage biased with a bias much larger than the temperature. 

Supplementary Discussion: Introduction to the basic concepts of transport 
statistics 

Central to the discussion of transport statistics is the concept of the probability distri- 
bution P(Q, At) of transporting Q charges in a time At. Defining the moment generating 
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function G of the distribution 

G(X,At) = e XQ P(Q, At), (7) 
Q 

the fcth moment (Q k )(At) of the distribution P(Q, At) is given by simply taking derivatives 
with respect to A, the generating variable, 

(Q k )(At) = d k x G(0) = £Q fc P(Q, Af). (8) 

Q 

If k — 0, then the result is 1 because the distribution is normalized. We further define the cu- 
mulant generating function as S{\ At) = logG(A, At), so the cumulants of the distribution 
((Q k )} are defined by 

((Q k ))(At)=d k S(0,At). (9) 

For example, «Q» = (Q),{{Q 2 )) = <(Q - (<3)) 2 ), ((<2 3 )> = <(Q " (Q)f), ({Q 4 )) = <(Q " 
(Q)) 4 ) — 3((Q — (Q)) 2 ), etc. The distribution may be expressed in terms of the cumulant 
generating function as 

P(Q, At) = J^- expf-AQ + 5(A, At)] (10) 

where the integration is over [0, 27ri]. 

For a random transport process to be Markovian means that the time evolution at time 
t + At is only determined by the state of the system at time i, where At is larger than 
any microscopic correlation time scale. Markovian noise is stationary if the probability of 
transmitting charge Qi in time Ati followed by charge Q 2 in time Ai 2 through a conductor is 
simply given by the product of the probability distributions, so the two events are statistically 
independent in time. This statistical independence implies that the cumulant generating 
function for transporting Q charges in time At = Ati + At 2 obeys the relation 5(A, At) = 
5(A, Ati) + £(A, At 2 ). This equation implies that the charge cumulant generator 5(A, At) 
must be linear in At if the Markovian random process is stationary, 

S(\,At) = AtH(X). (11) 

We define if (A) as the generating function of the current cumulants. The reason cumulants 
of current play a central role in the statistical characterization of transport, is the fact that 
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the cumulants of charge must also be linear in time. Therefore, the cumulants of current 
((I k )) are a time-independent characterization of the statistics of the transport, given by 

((I k )) = ((Q k ))/At = d k x H(0). (12) 



Supplementary discussion: Details of the Mathematical Derivations 

In order to describe the transport statistics of the interacting QD/QPC system, we begin 
with the simplest case of no transport through the QD. The statistical properties of the 
QPC are described by the current cumulant generating functions H a (X) = X k ((I k ))/k\^ 
where a = 1,2 indicates whether the QD is occupied or empty. These functions describe 
the intrinsic fluctuations of the QPC in the absence of any switching. 

To access the statistics of the QD and QPC when switching with rates Ti 2 is allowed, we 
now derive the joint generating function of both conductors, based on conditional master 
equation formalism. From the model described in the main text, transport through the 
quantum dot occurs through electrons tunneling from the left reservoir, onto the QD, and 
off through the right reservoir. If the number of electrons in the right (left) reservoir on 
either side of the QD is taken to be N R (N L ), then whether the current / is either Ii or I 2 is 
controlled by the resident charge of the QD, while the current through the QD is controlled 
by the transfered charge. The transfered charge M is given by counting the change in 
the symmetric combination (N R — N L )/2 from t = 0, while the resident charge is given 
by the change in —(N R + N L ) from t = (we note that it is important to choose these 
two independent variables to avoid introducing superfluous correlations into the system). 
We introduce the joint probability P a (N,M,t) of transmitting TV electrons through the 
QPC, and M electrons through the QD in the time t in state a. The mixed generator 
U a (\, M, t) = J2n eXN Pa{N, M, t) satisfies the master equation 



This equation describes the time rate of change of the probability of occupying state a in 
terms of incoming and outgoing contributions. The current transfered through QPC state 1 



ftWi(M) = (/Ti-ri)Wi(M) + r 2 W 2 (M + l/2), 



(13) 



d t U 2 {M) = (H 2 -T 2 )U 2 (M)+T 1 U 1 (M + 1/2). 



(14) 
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is counted by the first term in (13), while the possibility of switching to state 2 is accounted 
for in the second term (13) where N L decreases by one, accompanied by a change in M. 
Analogously, the current transfered through QPC state 2 is counted by the first term in 
(14), while the possibility of switching to state 1 is accounted for in the second term (14) 
where N R increases by 1, accompanied by a change in M. 

Introducing the full generator U a {\ x> t) = Em eXM ^a(X, M, i), the full conditional mas- 
ter equation is given by the matrix equation d t U = H[7, where 

y r ie */ 2 H 2 {\)-Y 2 ) 

For times larger than T^, the systems relaxes to the stationary distribution, controlled by 
the largest eigenvalue i7(A,x), so that ^2 a U a ~ exp[tH(\,x)]- Considering just the noise 
power (second cumulant) of the QPC alone, d\H{Q, 0), the result is 

«/ 2 » = E « /2 »« p « + (A^) 2 rir 2 /r^, (ie) 

a=l,2 

where Pi^ = r 25 i/(2r ) are the probabilities of occupying state a ) and ((I 2 )) a = H£(0) is 
the microscopic (shot) noise of the QPC in state a. The comparatively long switching times 
implies that the telegraph contribution to (16) dominates over the QPC shot noise, and 
therefore we may neglect the QPC shot noise. This amounts to replacing H a (X) —> XI a) 
which recovers Eq. (4) of the main text. 

From the definition of the generating function, the probability distribution of transmitted 
charge may be extracted by Fourier transforming on both variables, 

P(N,M,t) = J ^exp[tH(X, X )-xM-XN], (17) 

where the integration is over the interval [0, 27rz]. Replacing N — It,M — Jt gives the 
distributions of the two currents. In order to consider the stationary limit, it is necessary 
to have t > T^, and we may therefore evaluate the integral (17) in the stationary phase 
approximation. The dominant contribution to the joint probability distribution is then given 

by 

logP(7, J) = tmin A;X [tf (A, x) - JA - J X ], (18) 
which takes the form of a Legendre transform, yielding the result 



\ogP/t=-— (7-Jo)-r -2J[log(2J/fi)-l], fi = X_L2 > /(j_/ 1 )(/ 2 _/). (19) 
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The term inside the logarithm is responsible for the correlation between the two currents. 
The mixed generating functions may be calculated by only Fourier transforming on one of 
the above variables. Integrating over the A variable only gives the first mixed form, 

Hx{I,x) = - (gl 2 A/ 2)2 +fl ( eX/2 " ^ W 

where <? lj2 = y/Ti^I ~h~i\- The result in the main text Eq. (5) is given by normalizing the 
conditional distribution by subtracting off the X = contribution. Equation (20) recovers 
the result Eq. (3) in the main text for the special case of x = 0. The other conditional 
generating function is found by transforming on \ on ly> yielding the second mixed form 

w 2 (A, J) = xio - r + s - jio g [2 js/rir 2 ], (21) 

where S = J + ^/ J 2 + (AAJ — AT) 2 . The conditional cumulants (6) of the main text are 
found by taking derivatives of this generating function with respect to A. 

Applying the remaining Fourier transformation to either of the above forms yields the 
result (19), completing the circle. We note that the conditional distribution functions may 
be obtained by Fourier transforming with respect to either Hi(I,x) ~ 0) to obtain 

P(J|7), or with 7^ 2 (A, J) — 7^2(0, J) to obtain P(/| J), as a consequence of Bayes' theorem. 
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